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ABSTRACT 

The extreme towline tension prediction methods used by the U.S. Navy are 
based on analytical models developed by Milgram, Triantafyllou, Frimm and 
Anagnostou (1988) . The analysis of the test results for an instrumented offshore tow 
conducted in 1989 indicated that the analytical model underpredicted extreme tension. 
The reason for the underprediction of the extreme tension was the analytical model 
did not account for the second order ship surge motions and resulting low frequency 
tension fluctuations. Frimm and Milgram (1991) proposed a semi-empirical 
correction to the extreme tension model based on an analysis of the at sea test results. 

This thesis is a continuation of the above efforts to improve the prediction of 
extreme towline tensions in offshore towing. The goals of this thesis were to: 

1. Design and conduct a test program on the tow of a defueled nuclear submarine the 
ex-USS Ray by the salvage ship USS Bolster to determine the following: 

a The ability to predict the low frequency towline tensions; 

b. The ability of an in-line synthetic spring to reduce dynamic tensions; and, 

c. The accuracy of the semi-empirical correction to the predicted extreme 

tension. 

2. Develop and validate a numerical time simulation of the towing model from which 
the extreme tension statistics can be determined directly. 

The towing test program was successfully performed. The presence and 
significance of the low frequency tension variations were confirmed and 
improvements to the semi-empirical correction were developed. A numerical time 
simulation of the towing model was developed and tension statistics based on the 
simulated tension time history were calculated and compared with statistics generated 
from the test results and the existing analysis methods. The ability of an in-line 
synthetic spring to reduce dynamic tension was demonstrated. 

Thesis Advisor: Jerome H. Milgram 

Title: Professor of Ocean Engineering 
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Chapter One 
Introduction 



1.1 Motivation 

Modem seakeeping theory and the principles of extreme statistics can be used 
to improve the safety and efficiency of towing at sea (Milgram, Triantafyllou, Frimm, 
and Anagnostou, 1988). The recent loss of two oil barges off the coast of Puerto 
Rico due to towline failure is an example of the severity of the consequences of 
towline failure. The ability to accurately predict towline tensions under various 
environmental conditions influences towline design and the selection of towing speeds 
and courses (Milgram, 1993). 

The United States Navy is presently decommissioning nuclear submarines at a 
rate of approximately five submarines a year. When a nuclear submarine is taken 
out of service (decommissioned) she is towed, on the surface, from her homeport to 
Bremerton, Washington and scrapped. Though her nuclear fuel is removed, the 
reactor is in place, and a defueled nuclear submarine is still a high value asset. 

For many years the only approved method for towing U.S. Navy submarines 
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has been with an all wire towline attached to an Automatic Towing Machine (ATM). 
An ATM is mounted on the stem of the towing vessel and is designed to maintain 
constant tension in the towline. The ATM maintains a preset value of tension in the 
towline by paying out the towline when the tension exceeds the preset value and 
retrieving it when the tension is low. The number of United States Navy vessels 
equipped with ATM’s is decreasing at a rapid rate. In 1994 there will not be enough 
vessels with ATM’s to support the number of planned submarine decommissionings. 



In order to tow the defueled submarines from the shipyards where their fuel is 
removed to the location where they will be scrapped the U.S. Navy must have an 
alternative to all wire towlines and ATM towing. 



A typical towline consists of a short length of steel chain in series with a 
length of wire towline (Figure 1.1). 
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A towline’s stiffness is a function of both the towline geometry and the 
towline’s material properties. Endpoint dynamic motions introduce towline-water 
interactions which affect the relative influence of geometry and material properties on 
the total towline stiffness. 

A "stretcher" can be added to the towline as an alternative to the ATM for 
controlling towline tension. A stretcher is a length of large diameter fiber rope with 
low stiffness. As its name implies, a stretcher accommodates changes in endpoint 
motions by stretching. The low stiffness minimizes the change in towline tension for 
a given change in endpoint separation. A towline with a length of stretcher, wire 
towline, and anchor chain is called a composite towline. Figure 1 .2 shows an 
example of a typical composite towline. 



Composite Towline Configuration 



Towed Vessel 



. 90 feet of . 



Steel Chain 



_ 300 - 600 feet of . 

Synthetic Line 
(14 inch circumference) 







1000- 1500 feet 
— of Steel Wire — 
(2 1/4 inch diameter) 






> 



Figure 1.2 
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ATM’s are widely used in the U.S. Navy. Stretchers are common in the 



commercial towing industry. The purpose of this research is to evaluate the 
suitability of composite towlines for U.S. Navy use and avoid backfitting of ATM’s 
on the remainder of Navy’s fleet of ocean-going tugs. The suitability of composite 
towlines for U.S. Navy is to be demonstrated by accurately predicting the statistics of 
their extreme tensions in a range of configurations, applications, and environmental 
conditions. 

1 .2 Background 

The extreme tension is defined by the U.S. Navy Towing Manual (1988) as 
the sum of the steady and dynamic tension which has one chance in one thousand of 
being exceeded in one day of towing. The generation of the statistics of the extreme 
tensions in composite towlines is a mathematically and physically complicated 
problem due to the extremely non-linear behavior of the towline and the weakly non- 
linear hydrodynamic ship motions (Milgram, 1993). The behavior of this coupled 
non-linear system is driven by deterministic excitations- engine, throttle, and rudder 
changes) and stochastic excitations- the wind and the waves (Milgram, Triantafyllou, 
Frimm, and Anagnostou, 1988). The stochastic excitation of primary importance in 
towing is the excitation due to the wave-ship interactions. 
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Present towline extreme tension predictions use linear sea-keeping theory to 
predict ship motions. In linear sea-keeping theory the body’s response in each of its 
degrees of freedom to excitation is expressed as a frequency dependent linear transfer 
function (Newman, 1977). When a random process is operated upon by a linear 
system, some aspects of the statistics of the random process are preserved. Wave 
excitation is typically modelled as a gaussian process, and in linear ship motion theory 
the ship responses are gaussian. 

For most ship-sea interactions, linear sea-keeping theory works well. Non- 
linear wave-body interactions, specifically, those occurring at the difference 
frequencies of the incident irregular, random wave field, are important in open ocean 
towing (Faltinsen, 1990). Unlike motion in the pitch, heave, and roll modes of a 
freely floating object, surge motion does not experience large hydrostatic restoring 
forces. Therefore the surge motion of a vessel can be quite large. Normally, this is 
of little concern. Most sea-keeping prediction programs (MIT Five Degree of 
Freedom Sea-Keeping Program, 1975, and the David Taylor Model Basin Ship 
Motion Program, 1981), for example) do not predict surge motions. However, for 
two vessels coupled through a towline, large surge motions can cause large towline 
tensions. Accurate prediction of both first and second-order ship responses in 
irregular, random seas is important for this application and is one focus of current 
ocean engineering research. 
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The towing vessel and the towed vessel are free to move in all six of their 
respective degrees of freedom. The two vessels are coupled through the towline. The 
geometry of the system is shown in Figure 1.3. 



Towing Geometry 



Heave 




Figure 1.3 
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The ships’ motions are driven by the random wave excitation and time-varying 
cable tension. The full problem is a twelve-degree-of-freedom system with non-linear 
feedback (Milgram, et. al., 1988). 

The work of Milgram, Triantafyllou, Frimm, and Anagnostou (1988) and 
Frimm and Milgram (1991) form the basis for current methods for predicting the 
statistics of the extreme towline tensions in the offshore towing problem defined 
above. The former introduced the application of modem linear sea-keeping theory to 
the prediction of ship motions under tow. 

Frimm and Milgram (1991) used off-shore towing test data to supplement the 
previous work. The previous work was found to under-predict the extreme tensions. 
Frimm and Milgram (1991) used an empirical correction factor based on measured 
low-frequency tension data to compensate for the under-prediction of the extreme 
tensions. 

1.3 Approach 

The purposes of this research are to demonstrate that a composite cable with a 
stretcher is effective at reducing dynamic tension and to improve the accuracy of 
prediction of extreme towline tensions. 
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The effectiveness of composite cables with stretchers in reducing dynamic 
tension will be demonstrated through the evaluation of test data obtained during a 
major full-scale open ocean tow measurement program designed and conducted as 
part of this research. The extreme tension prediction methods of Milgram and Frimm 
will be compared with the test data. 

Extreme tension prediction methods will be improved by accomplishing the 
following: 

(1) Development of a numerical model for solving the cable dynamics equation for 
a composite towline; 

(2) Incorporating recent advancements in the prediction of second-order ship 
motions in the prediction of extreme tensions; 

(3) Performing full time simulations of the tug-tow system and determining 
extreme tension statistics from these time simulations; 

(4) Comparing the tensions predicted by incorporating the above improvements 
with the existing extreme tension prediction methods and the test data obtained 
during the open ocean tow; and, 

(5) Development and analysis of experimental apparatus for measuring relevant 
tow system parameters. 
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Chapter Two 
Theory 



In this chapter the underlying physical and statistical processes involved in the 
prediction of the extreme tension in a towline are presented. 

Problem Formulation 

The prediction of the extreme tension in a towline in open ocean towing requires an 
accurate characterization or estimation of the following: 

a. Sea state (frequency distribution of wave energy); 

b. Linear ship motion frequency responses; 

c. Second order ship surge frequency response functions; 

d. Towline tension changes due to dynamical ship motions; and, 

e. Towline extreme tension statistics using item d. 

The time varying excitations causing the ship motions are the wave-ship interactions 
and the ship-towline interactions. 
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Complications 



The numerical prediction of extreme tension is complicated by the following: 

a. The towline behavior is highly nonlinear and a solution of the cable dynamics 
equations can only be done numerically; 

b. The underlying physical process, the sea elevation spectrum is a random 
process; 

c. Prediction of second order ship responses is difficult and at the forefront of 
current research in ship hydrodynamics; and, 

d. The prediction of the response of a nonlinear system (the towline) operating on 
a random process (the sea elevation) is complicated and portions of the work 
must be performed numerically. 

2.1 Towline Dynamics 

The cable dynamics equations are (Triantafyllou, 1987): 

m T-7? = (T+T) (a + T-^) ~ b 4?4?~ T a (1) 

6t 2 5s 2 ot 5t 

with the dynamic tension, T, expressed as: 
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( 2 ) 



L 



L 



T = EA [ -^2 - — 
L L 




o 



0 



where: 



m = cable mass per unit length 
q = normal motions along cable 
T = static tension 
T = dynamic tension 
a = catenary static curvature 
s = Lagrangian coordinate along the cable 
b = sectional drag factor 
E = Young’s modulus 
A — cross-sectional area 

p 0 = tangential displacement due to mean static tension 
L = un stretched cable length 



These equations can not be solved analytically for the dynamic tension so they 
are solved numerically. The numerical solution is used to generate cable tension 
representations for use in determining the effect of the cable tension on the vessels in the tow 
system and for generating the extreme tension statistics. The two tension representations are 
a third-order polynomial representation and an equivalently linearized representation. 
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Polynomial Representation of Dynamic Tension 



For statistical calculations and for fast numerical towing simulations, it is convenient 
to represent the towline tension as a polynomial in terms of the towline extension and its 
time-derivative. Milgram, et. al. (1988) used a third-order polynomial. This representation 
is expressed as: 



T(S,e>=£ £ (3) 

m=Q n=Q 

with: 

m+nz 3 (4) 

and: 

(m, n) * (0, 0) (5) 

£(t) is the towline extension. 

The coefficients for the polynomial are found by minimizing the least square error 
between tensions calculated by the polynomial representation and tensions determined by 
numerically integrating the cable dynamics equations for a large set of cable extension time 
simulations (Milgram, et. al, 1988). The numerical integration scheme is a finite difference 
scheme with forward Euler time integration. 
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Figure 2.1 compares the results of the polynomial approximation with a measured 
towline tension. The measured towline extension, and its calculated time derivative, were 
used as inputs to the polynomial approximation. 




Figure 2.1 
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The polynomial approximation compares favorably with the measured tension data. 
The polynomial approximation generally over-predicts the peak tension. Thus, the 
polynomial’s use in predicting extreme tensions is conservative on the side of safety. 

In the experiment, the extension data was obtained by the use of a laser range finder. 
The tug-tow separation distance was measured during the data collection run and the 
difference between the measured separation and the mean separation distance during the run 
was recorded as the towline extension. 



The towline extension data used as input for the above analysis was numerically 
corrected to overcome measurement discontinuities which are described in detail in 
subsequent chapters. The recorded towline extension data is evenly spaced in time with a 
sampling rate for the laser range finder system of 2Hz. The laser range finder was hand- 
held on the towing vessel during a data collection run and focused on a reflective target on 
the submarine. During significant sea states the operator occasionally was unable to maintain 
the laser range finder "on target" for 1-2 second intervals. These "dropouts" complicated the 
data analysis. The data collection runs lasted 45 minutes. During a data collection run the 
operator would average 3 to 4 dropouts per second. The gaps in the data were filled in by 
a linear function between the endpoints of the dropout. These gaps affect the analysis since 
the elongation and its time derivative are used in the polynomial approximation. 
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Equivalently Linearized Cable Model 



Milgram et. al. (1988) used the method of equivalent linearization to develop a linear 
towline model for determining the influence of towline forces on the tug and tow motions. 
The method of equivalent linearization requires that the cross correlation function for the 
nonlinear system (RniXO) equal the cross correlation function for the equivalent linear system 
(R L (t)). That is: 

E[T NL (t) £(t + t) ] = E[T L (t) l (t + t) ] (6) 

where: 

Tnl is the nonlinear towline tension model in the form of the third-order 
polynomial 

T l is the equivalent linear towline tension model 
The equivalent linear towline tension model is to be expressed in the form: 

t l = * b„e it) (7) 

Milgram et. al. (1988) assume that the towline extension is a gaussian random process 
and show that: 

k eq = a 10 + 3a 30 m 05 + a 12 m 25 ( 8 ) 

and 
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a^j are the polynomial coefficients 

m 0? , m 2{ are the order zero and order two moments of the 
extension spectrum 

In this project, the equivalent linear towline model is used to characterize the 
relationship between towline forces and the vessels’ pitch and heave ship responses. This is 
permissible because the towline has only a minor influence on heave and pitch. The 
polynomial towline model is used in characterizing the effect of the towline force on the 
surge motions. 

Quasi-Static Representation of the Dynamic Tension 

Although the dynamic tension is most accurately determined by numerical integration 
of the cable equations or by use of the resulting polynomial approximation and the 
equivalently linearized model, a rough estimate and some understanding of the nature and 
importance of the cable dynamics can be obtained by using the cable’s static tension curve 
and measured extension time history. A plot of cable static tension as a function of extension 
can be generated from the solution of the equations of static equilibrium for a cable 
supported between two endpoints. Figure 2.2 shows a differential element of a cable and the 



forces acting on the element. 



T ♦ dT 




The equations of static equilibrium (Triantafyllou, 1990) are: 



-T - wsinOdl + (T + dT) cosd0 = 0 



(10) 



and 



-wcosOdl + (T + dT) sindO = 0 



(ll) 



with 



dx 2 + dz 2 = dl 2 



( 12 ) 



and 




( 14 ) 



The strain can be defined as: 



. dl - ds 
ds 



( 13 ) 
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For small angles, the above relations can be grouped to give the shallow sag cable 



equations (Triantafyllou, 1990): 



= wsinO (15) 

ds 



and 



~ = (1 + e) sin0 
ds 



(16) 



and 



T = WCOS0 

ds 



(17) 



and 



4^ = (1 + e) cos0 
ds 



(18) 



The cable equations can be combined and integrated (Triantafyllou, 1990) to yield: 



and 



T 



H 

COS0 



(19) 



tan0 = — - — (L - s) 
H H 



( 20 ) 



Figure 2.3 is a plot of the static tension-extension curves for a typical towline. 
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S tilde Tsnwon Exi*n*»on 




Figure 2.3 



This curve and the measured towline extension from a data run were used to generate 
the quasi-static tension versus time curve shown Figure 2.4. This quasi-static curve, when 
compared to the measured towline tension time history for the same extension time history, 
shows considerable estimation error. This demonstrates the importance of the cable 
dynamics in estimating towline tensions. 
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Quasi-Static - Dynamic Tension Comparison 




Time (seconds) 
Figure 2.4 



Towline Mechanical Properties 



Towlines used in U.S. Navy open ocean tows use wire, chain, and synthetic line in 
various combinations depending on the type of vessel to be towed. A summary of the 
important mechanical properties of each of these towline components are provided below. 

Wire Rope 
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Current U.S. Navy towing hawsers are usually 6 x 37 class, pre-formed, right-lay 
Improved Plow Steel (IPS), drawn galvanized wire rope of three types (U.S. Navy Towing 
Manual, 1988): 

a. 2-inch diameter, fiber core 

b. 2 1/4-inch diameter, fiber core 

c. 2 1/4-inch diameter, Independent Wire Rope Core (IWRC) 

Table 2-1 (U.S. Navy Towing Manual) summarizes the material properties of these 
wire ropes. 



Table 2-1. Wire Rope Characteristics 


Wire Type 


Diameter 


Weight in Air* 


Breaking Strength 




(inches) 


(lbs/ft) 


(lbs) 


Fiber Core 


2 


6.72 


288,000 


Fiber Core 


2 1/4 


8.51 


360,000 


IWRC 


2 1/4 


9.36 


387,000** 



* To get weight in water multiply by 0.87 

** Breaking strength for extra improved plow steel is 444,600 lbs 
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Fiber Rope 



Synthetic fiber rope acts as a spring in a composite towline. It responds to changes in 
dynamic loading by stretching with a relatively small change in tension. The U.S. Navy has 
used synthetic fiber rope in towlines for many years with mixed results. The synthetic fiber 
ropes originally used in the U.S. Navy for towing were made of nylon. Based on a large 
number of failures of nylon ropes, the U.S. Navy placed restrictions on their use (U.S. Navy 
Towing Manual, 1988) and initiated a program to evaluate the suitability of polyester towing 
hawsers. While nylon exhibits superior mechanical properties when dry, its properties when 
wet are significantly degraded. Table 2-2 (U.S. Navy Towing Manual, 1988) summarizes 
the relative characteristics of synthetic fiber ropes. 
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Table 2-2. Synthetic Fiber Ropes 


Fiber Type 


Strength 


Cyclic 


Bending 


Abrasion 


Heat 


Creep 






Fatigue 


Fatigue 


Resis. 


Resis. 






(1) 


(2) 


(2) 








Nylon (Dry) 


VG 


VG 


G 


E 


G 


G 


Nylon (Wet) 


G 


F 


F 


F 


- 


G 


Polyester 


VG 


VG 


VG 


VG 


G 


VG 


(Dry) 














Polyester 


VG 


VG 


G 


G 


- 


VG 


(Wet) 














Polypropy- 


F 


F 


P 


P 


P 


F 


lene (Dry) 














Polypropy- 


F 


F 


P 


F 


- 


F 


lene (Wet) 















E = Excellent VG = Very Good G = Good F = Fair P = Poor 

Table 2-3 summarizes the material properties of a sampling of the synthetic fiber 
towlines currently used by the U.S. Navy (U.S. Navy Towing Manual, 1988). 
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Table 2-3. Synthetic Fiber Rope Material Characteristics 


Fiber Type 


Weight in Air 


Circumference 


Tensile Strength 




(lbs/ft) 


(inches) 


(lbs) 


Polypropylene 

Three-Strand 

MIL-R-24049A 


2.33 


10 


115,000 


Dry Nylon 

Double-Braid 

MIL-R-24050C 


5.24 


14 


508,000* 


Polyester 

Double-Braid 

MIL-R-24677 


6.46 


14 


527,000 



*Wet nylon strength is 15% less than dry strength. 



Chain Pendant 



The chain pendant is used to connect the towline to the attachment point on the tow 
and serves one or more of the following functions (U.S. Navy Towing Manual, 1988): 

a. A chafing-resistant strong terminal connection to the towed vessel; 
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and, 



b. A means of absorbing dynamic loads in the towline, by virtue of its weight, 

which increases the depth of the catenary in the towline. 

Table 2-4 is a summary of the chain mechanical properties. 



Table 2-4. Anchor Chain Mechanical Properties 


Weight per Foot 


39.2 lbs/ft 


(in air)* 




Typical Length 


90 feet (one shot) 


Cross Sectional Area of Links 


6.28 in 2 


Young’s Modulus 


30 x 10 6 lbs/in 2 


Proof Test 


322,000 lbs 



* Weights in air must be multiplied by 0.87 to obtain weight in sea water 
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2.2 Ship Motions and Seakeeping 



Sway. Heave. Roll. Pitch, and Yaw 

The equations of motion for a single ship subject to wave and towline forces 
(Milgram et. al., 1988) are: 

|-co e 2 [M + A (co e ) ] + io e B (<o e ) + c|Xe i&>at 
= F (o> e ) e it>at + F H (<o e ) e i<Jat 

where: 

co e is the encounter frequency 
M is the mass matrix 
A is the added mass matrix 
B is the damping coefficient matrix 
C is the hydrostatic restoring force matrix 
X is the displacement vector 
F is the wave excitation force vector 
F h is the towline force vector 

In this project no change to the method for predicting the sway, heave, roll, 
pitch, and yaw wave-ship interactions or the towline influence on these five degrees 
of motion was made relative to the work of Milgram et. al. (1988). The wave-ship 
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interactions in these five degrees of freedom are determined by the MIT 5-D 
Seakeeping Program (1975). The MIT 5-D Seakeeping Program is based on the ship 
motion theory of Salvesen, Tuck, and Faltinsen (1970). This ship motion theory is a 
linear theory based on the following assumptions: 

a. The responses of a ship to irregular waves are a summation of the 
responses to regular waves of the contributing frequencies. 

b. Strip theory 1 accurately predicts ship responses to wave excitation. 

c. Viscous forces are negligible relative to inertial forces. 

d. Wave slopes are much less than one. 

Surge 



The equation of motion for a ship being towed or one towing another ship in 
the presence of waves can be written in the form of Newton’s second law: 



(m + m a (oi) ) = --i-pSC D (co) V ( t ) 2 + 



V(t) 



+ T(t) + F w (t) 



( 2 ) 



'Strip theory assumes that the variation in flow in a ship’s cross-sectional plane is much 
greater than the variation in flow in the longitudinal direction (Faltinsen, 1990). the ship is 
divided into a manageable number of sections and the hydrodynamic coefficients (added mass 
and damping) are calculated for each section and combined to determine the total ship added 
mass and damping coefficients. 
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where 



m = mass of the ship 

m, = added mass; the added mass in surge is small in comparison to 
the mass, typically 3 - 7% of m, and relatively insensitive to 
frequency of oscillation. 

V = velocity 

p = density of water 

S = wetted surface area 

C D = drag coefficient 

P = propulsive power 

T = towline force 

F w = force due to incident waves 

-The first term on the right hand side of 
equation (1) is the fluid drag force. 

-P/V(t) is the thrust due to the propeller. 

We multiply (1) by V(t): 

v ( t) (m + m a ) ^ = --|pSC D V(t ) 3 + P + V(t)T(t) + V(t)F w (t) (3) 

Then the surge is expressed as a small perturbation about a steady mean 
velocity: 
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V ( t ) = v c + € ( t ) 



( 4 ) 



We then rewrite (3) as: 



(V Q + e ( t) ) (m +m a ) d < v ° * e(t) \ = -ipSC 0 (V o + e(t) ) 3 + 

UL Z (5) 

p + (V c + e(t) )T(t) + (V 0 + e(t) ) F w (t) 



Expanding (5), dividing through by V 0 , and keeping terms up to first order in e gives: 



(m + m a ) = -lpSC D (V 0 2 + 3e(t)v 0 ) + 
_L + T ( t ) + € (t)T (t) + <V° ~ l ~_ 6 ^ t ^ .j.F„(t) 



( 6 ) 



V„ 



V, 



V„ 



Note that e(t)/V < < 1 and (V 0 + e(t))/V 0 « 1. Therefore we can rewrite (6) as: 



^£!(m +m a ) = -ApSC D (V 2 + 3e(t)V c ) + 
+ T(t) + F w (t) 

* 



( 7 ) 



P/V c is the steady thrust from the propeller. l/2pSC D V 0 2 is the steady drag. 
For V(t) = V„ 4- e(t) these terms are constants of equal magnitude. Under this 
assumption, the surge equation becomes 
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( 8 ) 



d€ d (t) (m + m a ) = T(t) + F w (t) - ^pSC D 3e(t) V 0 



Wave Excitation 



The wave excitation, F w (t), consists of non-linear and linear wave-ship 
interactions. 

Linear Wave Forces 



The linear, or first order, surge exciting force is estimated by the Froude- 
Krylov force. The Froude-Krylov force is the force corresponding to the unsteady 
pressure induced by the undisturbed incident waves. The Froude-Krylov force 
ignores the effect of the ship on the incident waves. 

The Froude-Krylov force is calculated by integrating the wave pressure on the 
ship’s hull surface. The wave pressure is obtained by rearranging the linearized 
unsteady Bernoulli equation into the form: 
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( 9 ) 



P = 



AM 

p 6t 



where 

p = pressure 
p = fluid density 
4> = wave potential 



The numerical integration of (9) is performed using the method of Milgram, 
Triantafyllou, Frimm, and Anagnostou (1988). This method uses the Gauss theorem 
to transform the surface integral of the wave pressure over the hull surface to a 
volume integral of the wave pressure gradient over the ship’s immersed volume: 



f fk = // s (ip) ’ nds = ^v-|^ dv 



( 10 ) 



The Froude-Krylov force oscillates at the frequency of the incident waves. 



Second Order Wave Forces 



The non-linear wave exciting forces are evaluated to second-order. The 
solution of the second order problem results in force with a non-zero mean and forces 
at the difference and sum of the incident wave frequencies (Faltinsen, 1990). 
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In linear wave theory, the quadratic velocity term in Bernoulli’s equation is 
neglected. This term is retained in the derivation of the second order wave forces in 
irregular waves. The pressure due to incident waves can be expressed as: 

p = ~P (|£) + (V<J>) • (V*) (11) 

where: 

p = incident wave pressure 

p = fluid density 

<j) = velocity potential 

For purposes of developing the form of the second order forces, the incident 
wave velocity potential will be expressed as the sum of two potentials of slightly 
different frequencies: 

<J) = e klZ sin(k 1 x - o^t) + -5^2 e k2Z sin(k 2 x - co 2 t) (12) 

U 1 <"> 2 

where 

g = gravitational constant 
A = wave amplitude 
o) = wave frequency 
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k = wave number = u) 2 /g 



The wave velocity potential consists of three separate potentials- the incident 
potential, the diffraction potential, and the radiated potential. For the purpose of 
illustrating the origin and nature of the second order forces, the velocity will be 
expressed as the gradient of a single potential: 

V = V(J) = o) 1 A 1 e klZ cos (a^t + 5 X ) + o) 2 A 2 e k2Z cos (o) 2 t + 5 2 ) 



where 

5j = phase angle in x. 

Substituting the square of (13) into (11) results in the following expression for 
the pressure: 
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P ■ p -tl " 2 [ “> 2 ¥ * “* 



A 2 
2 ^2 



+ oj 1 2 a 1 2 cos (2co 1 t + 25 x ) + 



( 14 ) 



o) 2 2 A 2 2 cos (2<o 2 t + 26 2 ) + (o 1 « 2 A 1 A 2 cos [ («! - <o 2 ) t +(6 X - 5 2 ) ] + 
a> 1 o) 2 A 1 A 2 cos [ (o> 1 + co 2 ) t + (b 1 + 8 2 ) ] ] 



This second order expression contains terms not found in the linear theory: 

a. Steady interactions, proportional to the oj i 2 A i 2 /2 terms. 

b. Unsteady interactions, oscillating at: 

1. Twice the wave frequencies 

2. The sum of two wave frequencies 

3. The difference of two wave frequencies 



The interactions at the difference frequencies are important to towing. These 
interactions, for waves with slightly different frequencies, produce slowly varying 
("drift") forces and moments which may cause resonance in surge (Faltinsen, 1990). 
If the natural frequency in surge is close to a difference frequency, large amplitude 
surge motion will occur, generating large towline tension. 
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The second-order wave-ship interactions at the difference frequencies are 
called the slow-drift forces and moments. Various methods have been formulated for 
estimating the slowly varying wave-drift interactions. The method chosen for this 
project is outlined below. 

Second Order Force Estimation 



The first step in determining the slowly varying surge force is to calculate the 
second order operators, R(ojJ, for the range of ship-wave encounter frequencies of 
interest. R(oi e ) is the ratio of second order force to wave amplitude squared. 

The encounter frequency, , is defined (Newman, 1977) as: 

co = u 0 - U cosP ( 1S ) 

g 

where: 

oi 0 = wave frequency 

U = ship speed 

g = gravitational constant 

j8 = angle between wave and ship headings 

The second order operators used for determining the slowly varying surge 
forces were generated from the seakeeping program NIREUS (Sclavounous, 1991). 
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NIREUS determines the second order operators (added resistance operators) by using 
strip theory and calculating, for each two dimensional section of a hull, the second 
order pressure due to the incident, diffracted, and radiated velocity potentials. 

The second order wave forces are significantly affected by the relative vertical 
motion between the ship and the incident wave surface (Gerritsma and Beukelman, 
1972). NIREUS takes the relative vertical motion into account in its calculation of 
the second order operators. A large number of closely spaced (in frequency) second 
order operators are generated by curve fitting through the operators calculated by 
NIREUS. 



Figure 2.5 is a typical wave elevation spectrum (Newman, 1977). The lowest 
frequency in a standard wave elevation spectrum is typically 0.2 - 0.6 hz and is a 
function of sea state. The second order forces of practical importance are below this 
frequency. At frequencies above this cutoff, the second order forces, occurring at the 
difference frequencies are small compared to the forces occurring at the wave 
frequencies. 
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(Jj — 



Figure 2.5 

There is a second order operator and a wave elevation spectral density 
associated with each discrete frequency. The second order wave force time history is 
expressed as the following: 



F„(t) - f f R( 

]=1 k=l ^ 



( 16 ) 



where: 



F sv is the second order force 



, to., + OK . 

R(^L- *) 



is the average of the second order operators 
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for the two frequencies 



If Wj - <jj k is small, the second order force oscillates with a frequency below 
the wave frequencies shown in Figure 2.5. If a>j - is large, the force oscillates at a 
frequency in the wave frequency range. The magnitude of the second order forces in 
the wave frequency range is small relative to the first order forces, on the order of 10 
- 100 times smaller (Faltinsen, 1990). 

Figure 2.6 shows two waves of equal amplitude and slightly different 
frequencies. 




Figure 2.6 
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Associated with each of the two waves is a second order operator (R(gj c )). 
Figure 2.7 shows the wave which results from the interaction of the waves shown in 
Figure 2.6. 




Figure 2.7 

The envelope enclosing the wave in figure 2.7 varies with a frequency which 
corresponds to the difference frequency of the two waves. The amplitude of the 
slowly varying envelope amplitude is equal to the amplitude of the original waves. 

The slowly varying force, F sv , is obtained by applying the average of the 
second order operators for the two waves to the amplitude of the slowly varying 
wave. 
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Comparison of Wave Force and Towline Force Magnitudes 



Although it is important to include the towline force in calculating ship 
motions, since it is the only spring-like surge restoring force, it is helpful to 
understand that the general magnitude of the time-dependent part of the towline force 
is small in comparison to wave force magnitudes. 

Figure 2.8 is a plot of surge acceleration versus time for a typical data run . 



Run 2 Soq* Acc^r*l*on 




500 1000 1500 2000 



TVrw (Mcond*) 

Figure 2.8 

The maximum surge acceleration of the towing vessel (tug) for the above run 
is 2.50 ft/sec 2 . The mass of the tug is 133,000 lbm. The total force on the tug is 
approximately 333,000 lbf. The maximum total towline tension on this data run was 
64,000 lbf and the maximum dynamic tension was 30,000 lbf. The maximum wave 
force is about five times larger than the maximum towline force or about ten times 
larger than the maximum dynamic towline force. 
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A similar conclusion regarding the relative magnitude of the wave and towline 
forces was reached by Hara and Yamakawa (1994) by evaluating the distribution of 
the peak ship motions. In linear ship motion theory the peak ship motions follow a 
Rayleigh distribution. If the magnitude of the towline forces is small, the peak 
motions should follow a Rayleigh distribution despite their presence. Hara and 
Yamakawa produced experimental results which demonstrate the minimal effect of 
towline forces on ship motions. However, low frequency surge resonance must be 
considered as an exception. 

System Natural Frequencies 

The tow-tug system can be pictured as the spring-mass-dashpot system shown 
in Figure 2.9. 




Figure 2.9 
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This system can be expressed as in matrix form as: 



M 1 0 


!> * 


bj o 


+ 


K 


-K 


0 m 2 


x 2 


0 b 2 


X 2 


-K 


K 



F x (t) 

F 2 (t) 



( 18 ) 



with: 



= vessel mass 
bj = vessel damping 
= towline stiffness 
F^t) = external forcing 

The system natural frequency is: 



, x 2 _ K M x +K M 2 

CO XT — 

M, M 2 



( 19 ) 



System Damping 



System damping is due to skin friction, wave radiation, eddy-making, and 
wave-drift damping (Faltinsen, 1990). 
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The system damping at wave frequencies (b;) is determined by the slope of the 



speed-resistance curve, 



b. = d (resistance) 
1 d (velocity) 



( 20 ) 



and, in the case of the tug, by the slope of the speed-resistance curve plus the 
propeller damping. Figures 2.10 and 2.1 1 are the speed- resistance curves for the tug 
and tow respectively. 



Resistance Curve for Tow 




Figure 2.10 
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Resistance Curve for ARS 38 




Figure 2.11 



The propeller damping is more difficult to calculate. The towing vessels used 
by the U.S. Navy are powered by diesel engines. For a given throttle setting , a 
diesel engine is a constant torque machine. As the vessel surges, the propeller inflow 
velocity changes, changing the effective thrust of the propeller. The propeller 
damping coefficient as a result of this motion is expressed as (Frimm ,1987): 



6T 

6v 



dK t _ dK 

2 p D 2 ( 1 - co ) 2 ^ [ dJ Q 

U 



—3 K t 
dJ 



2K q - 



J-^3 

dJ 



( 21 ) 
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with: 



T = thrust 
V = velocity 
p = water density 
D = propeller diameter 
w = wake fraction 
J = advance coefficient 
K, = thrust coefficient 
K q = torque coefficient 

For low frequency ship motions, the wave radiation damping can be neglected 
and, for low frequency motions in higher sea states, the wave-drift damping is the 
dominant damping component (Faltinsen, 1990). 

Stiffness 

The spring constants for a towline are found by solving the cable static 
equations, equations 2.1-19 through 2.1- 20 for the change in tension due to a change 
in end point separation. These equations were solved numerically for the composite 
cable. Figure 2.12 shows the tension-extension relationship for the composite cables 
evaluated during the submarine towing project. 



52 



Static Composite Towline Tension-Extension 




Figure 2.13 below is a plot of the spring constant as function of the towline 
extension for a typical towline. 



Composite Towline Spnng Constant 




Figure 2.13 
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2.3 Statistics 



Extreme Tension 
Theory 

The statistic of concern for U.S. Navy open ocean towing is the towline 
tension that has a 0. 1 % (a) chance of being exceeded in 24 hours (At) of towing 
(U.S. Navy Towing Manual, 1988). This tension is called the extreme tension. We 
denote this tension as: 

T a ,*t U> 

For a given time period the total tension consists of a steady component and a 
dynamic component. The steady component is governed by deterministic factors 
(towing speed, length of towline, tug and tow drag characteristics) while the dynamic 
component is a randomly distributed value. The statistical analysis required is the 
determination of the extremal statistics of the dynamic tension. 

For a given dynamic tension time history, the probability that a selected 
tension at time t+dt will be exceeded is a function of the tension and the rate of 
change of tension at time t (Newland, 1984). Figure 2.14, from Newland, shows this 
graphically. 
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The above concept is applied to the frequency with which the dynamic tension 
exceeds the extreme value. From the above it is clear that the joint probability 
density function (pdf) of the tension and its time derivative, p r r =/(T,T), is required 
to determine the frequency with which a tension level exceeds the extreme tension 
value. If the joint pdf is known, the rate of tension level crossings can be found by 
(Newland, 1984; Milgram et. al. , 1988): 

QO 

N(T 0 ) = J fp Ti f d(T, T) (2) 

o 
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If the peaks are distributed such that there is one extreme peak for each 
dynamic tension cycle (each zero level up-crossing), then the probability that the 
tension level will not be exceeded by any one tension peak is: 



( ^ext ) 



1 



N( 0 ) 



( 3 ) 



Milgram et. al. (1988) relate (3) and (1) by the following: 

a. The time interval. At, is equal to the product of the number of cycles, 
n, and the period of oscillation. The probability that the extreme 
tension will be exceeded during At is equal to : 

[p T (r 0 )] n ( 4 ) 

where P t (T 0 ) is the probability that T 0 will be exceeded for a single 
cycle. 

b. T a At is defined as the tension level which has a probability of a of 
being exceed during At. Substituting, (1) into (4) gives: 

= 1 - « < 5 > 

c. (5) and (3) are combined to yield the following: 
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( 6 ) 



r« fAt = Pf 1 [(l - a) "^T2 Tc] 

The validity of the above method for estimating T„ At is a function of the 
accuracy of the following assumptions: 

a. . The extreme tension peaks are statistically independent; and, 

b. There is only one peak tension for each zero up-crossing. 

The accuracy of these assumptions is determined by: 

a. The bandwidth of the underlying random process (Newland, 1984); 

b. The sample time, At,; and, 

c. The magnitude of the extreme peak. 

The accuracy of these assumptions can be verified by an analysis of the 
experimental and simulated tension time histories. 

Numerical Implementation 

The extreme tension statistics will be estimated two ways. The first way will 
be to use directly the probability computations of Milgram ,et. al. (1988) with an 
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empirical correction for the slowly varying ship motion and tension (Frimm and 
Milgram, 1991). The second way will be to perform a 1000 day time simulation of 
an open ocean tow and determine the extreme tension by histograms of the generated 
tension time history. The extreme tensions predicted by both methods- for the same 
vessels, towline, and environmental conditions, will be compared. 

Probability Computations 

This project is an extension of the work started by Milgram et. al. (1988) and 
continued by Christensen (1989) and Frimm and Milgram (1991). The procedure for 
predicting extreme tensions is taken from this earlier work with the improvements in 
solving the cable dynamics equations for composite cables and improved low 
frequency tension predictions developed in this project incorporated into the 
procedure. A full explanation of the numerical method for predicting the extreme 
tension is provided by Milgram et. al. (1988). The significant aspects of the method 
are: 

a. The rate of tension level crossings (2) is expressed as an expectation in 
terms of the pdf of the towline extension and its first and second time 
derivatives. This expression relies on the approximation of the tension 
as a third order polynomial in terms of the extension and its time- 
derivative (equation 2.1-3). 

b. The towline extension is approximated as a gaussian random process 
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with zero mean whose joint pdf, p { (£,£,£)> is expressed in terms of its 
spectral moments using, from (Parzen, i960),: 



U x , 



(2lt) n/2 ! K \ 1/2 



exp[-l £ (xj) Ki k (x k ) ] 

J t k = 



(7) 



where: 

K is the covariance matrix. The matrix elements are determined by the 
solution, in the frequency domain, of the twelve degree of freedom 
model of the tug-tow system. 

The accuracy of (a) is driven by the accuracy of the polynomial expression for 
the tension. 

The validity of (b) is determined by the extent to which the extension is a 
gaussian random process, which, in turn, is determined by the extent to which the 
ship motions can be considered a linear time invariant process operating on a gaussian 
sea state. 



Milgram, et. al. (1988) use strictly linear ship motion theory, justifying the 
use of (b). Second order surge motions, not included in the earlier theory have been 
deemed to be important (Frimm and Milgram, 1991). Since these second order 
motions lead to a non-gaussian process, including them requires an altered theory 
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(Longuet-Higgins, 1963). Frimm and Milgram (1991) concluded that a conservative 
estimate of the increase in total extreme tension due to the second order motions is 
twice the RMS value of the low frequency tension. 

Frimm and Milgram (1991) estimated the low frequency tension contribution 
by determining the low frequency tension component of a measured tension time 
history as a function of sea state. The measured data were extrapolated to higher sea 
states. Figure 2.15 is a plot of the estimated low frequency contribution to towline 
tension. 




Figure 2.15 
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Time Simulations 



The concept behind generating statistics from time simulations is 
straightforward. The twelve degree of freedom tug-tow system is modelled 
completely, solved in the time domain for the time period of interest, and the 
generated towline tension time history is scanned to determine the extreme tension for 
the time period. 

Like the method above, accuracy of the time simulations depends on an 
accurate representation of the wave-ship interactions, an accurate representation of the 
waves, and an accurate towline model. However, the accuracy of the time 
simulations relative to the method outlined above should improve due to the 
elimination of the probabilistic approximations. 

In the time simulations the vessels’ relative position, velocity, and acceleration 
are determined at time t Q by solving the vessel equations of motion with the excitation 
provided by the waves and the towline. The wave spectrum used in the simulations is 
either a measured sea spectrum from a data run or a Pierson-Moskowitz spectrum. 

The towline tension model used the equivalent linear model for linear ship heave and 
pitch motions and the polynomial tension expression for the surge simulation and 
extreme tension prediction. The vessel displacements are calculated form the original 
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position are calculated and used to update the polynomial and equivalent linear 



towline models. 



2.4 Towline Extreme Tension Prediction 



Problem Solution 

Two methods for predicting extreme towline tensions in offshore towline 
presently exist: 

a. Nonlinear extreme tension prediction program. 

b. Time domain simulation of the complete towing problem. 

Figure 2.16 shows the basic flow path for the nonlinear extreme tension 
prediction program (Milgram, et. al, 1988). Figure 2.17 shows the flow path for the 
time simulation. 
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Nonlinear Extreme Tension Prediction Program 
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Figure 2.16 
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Time Simulation 
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Figure 2.17 
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The computer programs required to complete the process charted in 
Figure 2.16 are: 

a. MIT 5D Seakeeping Program. MIT 5D requires as input: 

1. Tug underwater hull geometry, naval 
architectural characteristics, and propeller 
characteristics. 

2. Tow underwater hull characteristics, naval 
architectural characteristics, and propeller 
characteristics. 

3. Towing speed and relative ship-wave heading. 

4. Frequencies for calculating ship motion 
responses. 

NEWSRG1, a linear surge response program. NEWSRG1 requires 
similar input to that required for executing the MIT 5D program. 
CABNEXP, a program for determining the coefficients for the 
polynomial representation of the towline tension. CABNEXP requires 
as input: 

1. Towline segments’ material and geometric 
properties. 

2. Anticipated mean tension. 



b. 



c. 
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d. PMSPEC, a program for generating a Pierson-Moscowitz wave 

elevation spectrum 1 . The required input for PMSPEC is wind speed. 

d. 12DNEW, a program which solves the coupled 12 degree of freedom 
problem defined earlier, using output from each of the above programs. 

e. NLEXTR5, a program for predicting extreme towline tension. 
NLEXTR5 requires as input: 

1. Towline extension spectral moments, m 0 , m 2 , 
and m 4 (Output from 12DNEW). 

2. Mean tension, towline length, tow speed, 
relative ship- wave heading, and wind speed. 

3. Coefficients for the polynomial 
representation of the towline tension. 

The semi-empirical correction to the towline tension due to low frequency 
motions identified earlier is added to the extreme tension predicted by NLEXTR5. 

The programs required to complete the process indentified in Figure 2.17 are: 

a. MIT 5D Seakeeping Program. 

b. NEWSRG1. 



‘Any sea spectrum can be used as this portion of the input. A measured sea spectrum is 
preferred for some applications. 
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c. 


CABNEXP. 


d. 


PMSPEC. 


d. 


12DNEW. 



e. RWAVETOW, a program for determining the first and second order 
forces acting on the tug and tow vessels. RWAVETOW requires the 
following as input: 

1. Output from NEWSRG1 and 12DNEW. These 
outputs are the ship motion transfer functions. 

2. Second order surge added resistance operators 
(transfer functions). 

f. SURGESIM, a program for calculating the towline tension time history 
using the first and second order forces calculated by RWAVETOW. In 
addition to the output from RWAVETOW, SURGESIM requires as 
input: 

1. Tug and tow masses and added masses. 

2. Tug and tow damping coefficients. 

3. The coefficients for the polynomial 
representation of the towline tension. 

SURGESIM is run for 2500 days and the highest tension peak is read directly 
from the time record. 
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Chapter Three 



Experimental Preparation 



3.1 Background 

Introduction 

The data for this project were acquired during a "Tow of Opportunity". This 
was the Navy tow of the SSN 637 class submarine ex-USS Ray (SSN 653), first by 
the ARS 50 class salvage ship USS Grapple (ARS 53), and then by the ARS 38 class 
rescue ship USS Bolster (ARS 38). The tow of the ex-USS Ray occurred during the 
summer of 1993. Figure 3.1 is the track of the tow of the ex-USS Ray. It left 
Charleston, South Carolina on July 18, and arrived in Bremerton, Washington on 
September 5. 



69 




Puget Sound 



San Diego 



Charleston, S.C. 



Panama Canal 



Figure 3. 1 



The tow consisted of three distinct legs- from Charleston, South Carolina to 
Panama, from Panama to San Diego, California, and from San Diego, California to 
Bremerton, Washington. In anticipation of weak seas and anticipated equipment 
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debugging while in the Atlantic Ocean, the first leg was to be used as a calibration 



run. 



Since the tow was a scheduled Navy event the data acquisition had to be 
planned to have little or no impact on the progress of the tow. Additionally, 
personnel are prohibited from boarding nuclear submarines being towed. That is, 
once underway, the equipment onboard the submarine can not be adjusted or repaired 
until the tow is pierside. The above two factors significantly affected the data 
acquisition plan and the data acquisition equipment selection. 

Data Requirements 

A considerable amount of equipment is required to obtain the data necessary to 
fully quantify the twelve-degree of freedom problem outlined in the previous chapter. 
Data must be acquired in four areas. These areas are: 

1. Tug motions. 

2. Tow motions. 

3. Towline tension and elongation. 

4. Environmental (Wind and Waves). 
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Identification of Critical Parameters 



Prior to selecting the equipment required to support data acquisition in the four 
areas listed above, the parameters essential to the prediction of towline extreme 
tension statistics and to the validation of the analytical and numerical models must be 
identified. Based on a review of similar towing projects (Frimm and Milgram, 1991; 
Christensen, 1989), the following parameters were considered to be critical to the 
success of the project: 

1. Towline tension. 

2. Towline elongation. 

3. Wave spectra. 

Towline Tension 



The importance of accurately measuring towline tension is obvious. 



Towline Elongation 



Towline elongation is vital for two fundamentally important reasons. 



First, an accurate time history of towline elongation can be used to validate the 
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composite towline numerical model directly. This can be done in two ways. The 
first is to use the measured elongation as "input" to the cable differential equations 
and make a comparison between the tension time history predicted by the numerical 
solution and the measured tension time history. Alternatively, the towline elongation 
and its calculated time derivative can be used as the end point motion input to the 
numerical composite towline model and the predicted towline tension time history can 
be compared with the measured towline tension time history. 

The towline elongation can provide an estimate of the tow system’s behavior. 
For this project, the tow’s displacement and drag were significantly larger than that of 
the tug. The magnitude of the tug’s wave frequency surge motions was estimated to 
be significantly larger than that of the tow. Therefore, most of the towline elongation 
at wave frequencies was expected to be due to tug surge motion. The towline 
elongation provides a reliable measure of the relative surge motions between the two 
vessels. 

Wave Measurements 



Two principal options are available for measuring waves, ship-mounted 
sensors and wave measurement buoys. Ship-mounted wave sensors use a microwave 
radar system that bounces radar waves off the wave surface and measures the doppler 
shift of the returned signal (Christensen, 1989). The doppler shift is proportional to 
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the surface relative velocity. A fundamental problem with ship mounted wave 
measurement devices is that the ship and ship motions influence the surrounding wave 
field (Sellars, 1967). Additionally, the ship’s pitching and heaving motions must be 
known and the doppler shift corrected for these motions to get the actual sea motion. 

The wave motions in the proximity of a moving ship are the sum of the 
incident waves, the waves diffracted by the ship, and the radiated waves due to the 
ship’s motions. Ship mounted wave measuring devices are typically bow-mounted on 
extension arms to maximize the distance between the ship and the measuring device 
and thus minimize the ship’s influence on the measured waves. The wave 
measurement device is calibrated, and ship influence corrected for, by comparing the 
device’s output with the wave spectrum measured by a wave measurement buoy. 

Previous towline measurement projects suffered from the use of poorly 
constructed ship mounted wave measuring devices (Frimm and Milgram, 1991). The 
output of the ship mounted wave measurement device used by Frimm and Milgram 
was compared to the wave spectra measured by National Oceanic and Atmospheric 
(NOAA) buoys. Based on this comparison they determined that the spectra measured 
by the ship-mounted devices was unusable. 

Based on the experience of Frimm and Milgram wave buoys were used for 
wave measurement. 
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3.2 Equipment Selection 



A general description of the equipment required to obtain data in the four 
areas identified earlier follows. The selection of the equipment was strongly 
influenced by the work of Christensen, "Plans and Specifications For A Full Scale 
Towing Model Validation Experiment" (1989). A detailed description of the 
experimental apparatus is provided in Appendix A. 

Tug Motions 

The motions of the tug in its six degrees of freedom were measured using a 
"Six Degree of Freedom Stable Platform" leased from the National Research Council 
of Canada’s Institute for Marine Dynamics (IMD). The Six Degree of Freedom 
Stable Platform provides analog output proportional to the following: 

1. Roll 

2. Pitch 

3. Yaw 

4. Heave 

5. Surge 

6. Sway 

7. Roll Rate 
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8 . 



Pitch Rate 



9. Yaw Rate 



The above motions are defined in Figure 3.2. 



Heave 




Figure 3.2 
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It is desirable to locate the motion sensors near the center of gravity (CG) of 
the vessel. Installed shipboard equipment prevented locating the Six Degree of 
Freedom sensor at the vessels’ CG’s. The acquired motion data were corrected for 
the separation between the sensor’s installed location and the ship’s CG. 

The tug’s speed over ground was measured using a Global Positioning System 
(GPS) receiver. The GPS receiver provides a continuous record of the ship’s latitude 
and longitude. The accuracy of the receiver is +/- 30 meters. 

The tug’s heading partly defines the tow geometry. The tug’s heading was 
measured using an electronic fluxgate compass. 

Tow Motions 



A Six Degree of Freedom Stable Platform, a GPS receiver, and a fluxgate 
compass were used to acquire the motion, speed, and heading of the submarine. 
Acquisition of the tow data was made difficult by the fact that the tow could not be 
boarded after leaving port. Therefore, the tow had to be provided with an electrical 
power system and a system for transmitting data to the tug for analysis. 

The relative position of the tug and the tow, as measured by the two GPS 
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receivers, was to be compared with the tug heading to determine the time history of 
the "shearing" of the tow. When a tow deviates form a position directly abaft the tug 
it is said to be "shearing". The acquired GPS data were not accurate enough for this 
purpose. 

Electrical Power 



Table 3-1 summarizes the power requirements for the equipment onboard the 

tow. 



Table 3-1. Tow Vessel Power Requirements 


Item 


Current (Amps) 


Voltage (Volts) 


Power 








(Watts) 


Stable Platform 


2.5 


28.0 DC 


70.0 


Compass 


0.080 


12.0 DC 


0.1 


GPS Receiver 


0.2 


12.0 DC 


2.4 


FM Telemetry Modem 


4.0 


12.5 DC 


25.0 


Load Cell Signal 


0.001 


12.0 DC 


0.01 


Conditioner 
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Data Acquisition 


2.7 


120.0 AC 


324.0 


Computer 








Computer Monitor* 


0.7 


120.0 AC 


84.0 


Total 


9.98 


- 


505.5 



* Turned off during towing. 



Six marine gel cell batteries were used as the primary energy source for the 
equipment listed in Table 3-1. The cells provide 12 volt DC power and have a 450 
amp- hour rating. 

The batteries did not have sufficient stored energy for continuous operation of 
the data acquisition system and transmitter over the course of any of the three legs of 
the tow. Therefore, an automatic timer was installed to turn on the system only 
during the data acquisition runs. Data acquisition was limited to two 45 minute 
periods each day. 

The batteries were charged during towing by two solar powered and two wind 
powered electrical generators. The two solar powered generators were each rated for 
48 Watts @ 14 volts. The two solar powered generators were capable of producing a 
total of 36 Amp-hours per day. The two wind powered generators were each rated 
for 6 Amps at a wind speed of 19 knots. 
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The wiring diagram for the electrical system is shown as Figure 3.3. 




Figure 3.3 



Telemetry 

The digitized data acquired on the tow was transmitted to the tug via a Radio 
Frequency (RF) modem. The analog outputs from the data acquisition equipment 
were linked to the RF modem through a computer with an analog-to-digital converter 
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housed on the tow. 



The advantages and disadvantages of analog FM telemetry and digital data 
transmission via an RF modem are discussed fully in previous works (Christensen, 
1989; DeBord, Purl, Mlady, Wisch, & Zahn, 1987) and are summarized in Table 3- 
2 below. 



Table 3-2. Telemetry System Relative Merits 


System 


Advantages 


Disadvantages 


FM 


Transmits raw data- analysis not 
limited by frequency constraints 


Frequency sensitive 


Digital 


Tolerate frequency changes 


Bandwidth constraints 


Allows multiplexing 



The data acquired on the tow were also stored on a computer on the tow. 



Towline Tension and Extension 



Tension 
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The towline tension was measured by load cells installed 
at each end of the towline. The load cells are large links instrumented with strain 
gages for measuring towline tension. The specifications for the load cells are 
provided in Table 3-3. 



Table 3-3. Load Cell Specifications 


Item 


Specification 


Capacity 


200,000 lbs 


Sensing Element 


Dual Bridge 350 OHM Strain Gages 


Output 


2mV per volt excitation at rated capacity 


Excitation 


10 volts DC 


Calibration Error 


Less than 0.25% of full capacity 


Material 


Stainless Steel 



The installed configuration of the load cell at the tow end is shown in Figure 
3.4. 
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Figure 3.4 

The installation of the load cell at the tug end is more complicated. Normally 
the towline is at the tug end is attached to an automatic towing machine (ATM). The 
ATM controls the tension in the towline by paying out or retrieving towline as 
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necessary to maintain a constant towline tension. During the data acquisition period 
the towline is connected to the load cell using a carpenter stopper. Figure 3.5 shows a 
carpenter stopper (U.S. Navy Towing Manual). 




Figure 3.5 
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Figure 3.6 shows the installed configuration of the carpenter stopper and load cell. 




Figure 3.6 



A load cell was installed at each end of the towline for two reasons. 



Towline 
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tension is one of the three critical parameters for a successful test program and a 
second load cell provides redundancy. Additionally, concern was expressed by others 
over the possible existence of a delay between when a tension increase at one end of 
the towline is felt at the opposite end. 

Extension 



The importance of accurately measuring towline extension can not be 
overstated. The most accurate and reliable method for determining towline extension 
is to measure the tug and tow separation distance. The tug and tow separation 
distance provides the towline end point displacement time history. The end point 
displacement is then used in the cable equations of motion to determine cable 
extension. 

An infrared laser range finder was selected for measuring tug and tow 
separation. The advantage of the laser range finder over radar range finders is that 
the relatively short measurement intervals of the laser range finder are superior for 
unstable and moving platforms (Christensen, 1989). 

The laser range finder was hand-held. Laser-light reflective targets were 
placed on the submarine’s sail. The output of the laser range finder was fed directly 
to a data acquisition computer. 
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Continuous tug-tow range information was difficult to obtain in heavy seas. 
The operator often had difficulty "hitting" the target on the submarine (drop-out). 
When dropouts occur the data was smoothed to compensate. Figure 3.7 shows 
representative laser range finder output after corrective steps were taken to 
compensate for dropouts. 



SAMPLE LASER RANGE FINDER OUTPUT 




Figure 3.7 



Wave Measurement Buoys 



Background 

The problems encountered on previous projects with ship-mounted wave 
measurement devices and the operational constraints imposed by a "tow of 
opportunity" drove the selection of the wave measurement devices. Semi-disposable 
wave measurement buoys were designed, developed, and used for measuring the wave 
spectra. 

Theory 



An object floating on waves much longer than the object is subject to the 
dynamic conditions of the water it displaces (Principles of Naval Architecture, 1967). 
Froude ("On the Rolling of Ships", Institute of Naval Architecture Transactions, 

1861) confirmed experimentally that the sum of the forces acting on a particle on a 
wave’s surface is perpendicular to the wave’s surface by showing experimentally that 
a plumb bob suspended from a mast on a wave float would remain above the center 
of buoyancy of the float as the float rides the free surface (Figure 3.8). 
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Figure 3.8 
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Figure 3.9 is a free body diagram of the wave buoy. 




Figure 3.9 

The wave buoy was constructed using a single accelerometer, mounted such 
that the axis of the accelerometer is vertical when the buoy is horizontal. The 
accelerometer can only sense accelerations along its longitudinal axis. In order to 
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